Two-color stabilization of atomic hydrogen in circularly polarized laser fields 
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Dynamic stabilization of atomic hydrogen against ionization in high-frequency single- and two- 
color, circularly polarized laser pulses is observed by numerically solving the three-dimensional, 
time-dependent Schrodinger equation. The single-color case is revisited and numerically determined 
ionization rates are compared with both, exact and approximate high-frequency Floquet rates. The 
position of the peaks in the photoelectron spectra can be explained with the help of dressed initial 
states. In two-color laser fields of opposite circular polarization the stabilized probability density 
may be shaped in various ways. For laser frequencies Wi and ui2 = nui, n = 2, 3, . . . and sufficiently 
large excursion amplitudes n + 1 distinct probability density peaks are observed. This may be 
viewed as the generalization of the well-known "dichotomy" in linearly polarized laser fields, i.e, 
as "trichotomy," "quatrochotomy," "pentachotomy" etc. All those observed structures and their 
"hula- hoop" -like dynamics can be understood with the help of high-frequency Floquet theory and 
the two-color Kramers-Henneberger transformation. The shaping of the probability density in the 
stabilization regime can be realized without additional loss in the survival probability, as compared 
to the corresponding single-color results. 

PACS numbers: 32.80.Fb, 32.80.Rm, 42.50.Hz 



I. INTRODUCTION 



For sufficiently high laser frequency and intensity, the 
ionization probability of an atom decreases by further in- 
creasing the laser intensity. This, at first sight, counterin- 
tuitive effect is called dynamic or adiabatic stabilization, 
depending on the pulse length 0, | |, |, § § 0, |. 
So far, it has been observed in numerical simulations 
only (see || and, for very recent results, 0, |ll|). Its 
theoretical explanation is most elucidatingin the frame- 
work of high-frequency Floquet theory ||]: the time- 
averaged Kramers-Henneberger potential, as it is seen by 
a free electron moving in the laser field, provides bound 
states. The electronic probability density of these states 
is spread over spatial regions with dimensions of the order 
of the classical excursion amplitude which may be tens 
of atomic units or more. The laser frequency u>, required 
for dynamic stabilization where the electron starts from 
the field- free ground state with binding energy £q, has to 
exceed the field-dressed binding energy, hu> > \£q\. Since 
such high frequencies are not yet accessible to nowadays 
intense laser systems the experimental verification of dy- 
namic stabilization is still to come. However, with the 
intense, short-wavelength, and coherent light sources un- 
der construction worldwide (e.g., at DESY in Hamburg) 
ground state stabilization should become experimentally 
accessible soon. Signatures of dynamic stabilization of 



Rydberg atoms , where \£q\ is smaller, have been mea- 
sured |12| . 

Dynamic stabilization has been studied extensively 
during the last decade. Ionization rates and energy 
shifts were obtained numerically by Floquet theory |15| ] 
or approximate, analytical high-frequency Floquet the- 
ory j^, Q| . Full numerical ab initio solutions of the time- 
dependent, single-electron Schrodinger equation were ob- 
tained for low-dimensional model systems (see, e.g., 

@, n m, 0, oj, H) and ' after the early work in ® 

very recently, in three spatial dimensions for linear po- 
larization [jlO| and circularly polarized laser pulses [ pTj . 
Although the Floquet ionization rates decrease monoton- 
ically for laser intensity / — * oo the ionization probability 
in full ab initio solutions of the TDSE for finite laser 
pulses assumes a minimum and increases again there- 
after. In the limit of laser intensity / — > oo but other- 
wise fixed pulse parameters the survival probability of 
an atom will thus approach zero. Moreover, it has been 
found that the magnetic field of the laser — usually ne- 
glected because of the dipole approximation — leads to in- 
creased ionization [^0| . This is due to the ponderomotive 
force which pushes the electrons in propagation direction. 
In multi-electron atoms electron correlation also hinders 
stabilization (see J2l|, and references therein). Two-color 
interference effects on stabilization in linear polarization 
" "231. 



was studied in 22 
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In this paper the dynamic stabilization of H in cir- 
cularly polarized single- and two-color fields is investi- 
gated. The article is organized as follows: in Sec. II the 
three-dimensional, time-dependent Schrodinger equation 
to be solved numerically is introduced, and the form of 
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the applied two-color fields is explained. In Sec. Ill the 



results concerning the single-color, circular polarization 
case are presented. Besides snapshots of the probabil- 
ity density and photoelectron spectra, particular empha- 
sis is put on the comparison of ionization rates and en- 
ergy shifts with results from Floquet theory. In Sec. IV 
the main two-color results are presented. Probability 
density structures are observed which display well sep- 
arated maxima, the number of which is related to the 
frequency ratio of the two fields. The structure moves as 
a whole in the laser field and can be explained with the 
help of two-color, time-averaged Kramers- Henneberger 
potentials. Survival probabilities after the two-color 
pulses are presented as a function of the two laser's peak 
field strengths. The two-color survival probabilities and 
above-threshold ionization spectra are compared with the 
corresponding single-color results where the second field 
was absent. It is good news that adding a second color is 
possible without a drastic inhibition of stabilization even 
if the second laser's intensity lies in the "death valley" of 
high ionization. Finally, we conclude in Sec. [v|. 



is solved numerically in three spatial dimensions (atomic 
units are used unless noted otherwise) . The intense laser 
field can be treated classically. Since the laser wave- 
length is large compared to all other relevant length 
scales and the electron dynamics is non-relativistic for 
our laser parameters, the dipole approximation can be 
applied so that the vector potential A(t) does not de- 
pend on spatial coordinates. Hence, the magnetic com- 
ponent of the laser field is neglected. If we restrict our- 
selves to study coplanar two-color fields (of, however, ar- 
bitrary ellipticity) we may choose the coordinate system 
in such a way that the z-component of the electric field 
E{t) = -d t A(t) vanishes: E(t) = E^{t) + E^(t) 



(!)/ 



The wave function ^(r, t) is expanded in spherical har- 
monics, 



II. THEORY 

The time-dependent Schrodinger equation (TDSE) for 
the electron of atomic hydrogen in the electric field of a 
laser 
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i^*(r,i)= Q[-iV + A(i)] 2 -~ )*,,./! ,1, 



which leads to a set of coupled TDSEs for the radial wave 
functions (&e m (r,t), 



i^e m {r, t) 
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r 



is the effective atomic 



where Vf(r) = -± + « 
potential including the centrifugal barrier, A(t) = 
A x (t) + iAy(t), N± m = v / (eTm)(e±m + l)/2, and 
the purely time-dependent j4 2 -term has been trans- 
formed away. The matrix elements (£m\LM\£'m'} := 
/ d(cM&)d(pY£ m (>9,v)Y LM (>&,tp)Yi, m >('&,(p) maybe ex- 
pressed in terms of Clebsch-Gordan coefficients. 

It is worth noticing the coupling of different £ and m 
quantum numbers in the TDSE (||) due to the laser field: 
there is a £m — > £ ± 1, m ± 1 coupling as well as a £m — > 



£^1, m±l coupling (where compatible with the condition 
— £ < m < £). From the knowledge of these couplings one 
can immediately deduce non-perturbative selection rules, 
e.g., if the initial wave function is the Is groundstate, 
states with {£ = 1, m = 0), with {£ = 2, m = ±1), with 
{£ = 3, m = ±2) etc. are 'unreachable' and thus remain 
unoccupied. 

In the actual implementation, the radial coordinate r 
was discretized directly in position space. The details of 
the implementation may be found in |p3 and is similar 
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to what has been described in [B5| for the simpler case of 
linear polarization. 

The radius of our numerical grid, including the region 
where the non- vanishing imaginary potential damps away 
probability density approaching the grid boundary, was 
between 240 a. u. and 450 a. u. The maximum i quantum 
number was typically ^ max = 30. For testing the con- 
vergence of our results runs up to £ max = 50 were also 
performed — with no significant changes in the observ- 
ables of interest. For obtaining the results presented in 
this paper, the run times of our Schrodinger-solver were 
always less than 18 hours on 1 GHz PCs. 

We focus on two-color fields where both fields are circu- 
larly polarized. We assume a sin 2 -shaped up- and down- 



A/23 cycles an d 



ramp of the zth field, i = 1, 2 over iV{ 

a constant part of n[ 2 cycles duration. In the following 

we will refer to such a pulse as a (N^\ n[ 2 , )-pulse. 
For the ith electric field in, say, x-direction 



E. 



(i), 



■udW 

£?W cosojit 

« 



for < t < Ti, 
for T x < t < T 2 , 
E™(t-T 2 + Ti) for T 2 <t<T 3 , 
otherwise 

(3) 

is set where E^(t) = sin 2 (7rf/2Ti) coscjjt The ab- 
solute ramping time as well as the duration of the con- 
stant part of the laser pulse are chosen the same for 
both colors so that the times T\ = 2nN^ /<j,, T 2 = 

T\ + 2ttN^2 /u)i, etc. are independent of i. Since we are 
interested in opposite circular polarizations 

E?\t) = E^(t - 37T/2W0, Ef\t) = Ei 2 \t - ir/2^) 

(4) 

are taken as the y-components of the two colors. The 
corresponding vector potential can be calculated from 
A(t) — — J* dt' E(t'). In such two-color pulses, as long 

as — , are integer numbers, a free electron 
(initially at rest) returns to its starting point at the end 
of the laser pulse without residual drift momentum. 



III. THE CIRCULAR SINGLE-COLOR CASE 
REVISITED 

Before turning to the two-color study let us analyze the 
simpler single-color case. For the latter, results from ab 
initio solutions of the three-dimensional TDSE were ob- 
tained only recently [ |Toj] whereas two-dimensional model 
solutions were presented earlier @, || @. Ioniza- 
tion rates as well as energy shifts for hydrogen in the sta- 
bilization regime (with dipole approximation) are avail- 
able for several years |D| so that a comparison with our 
finite pulse results is possible and will be presented in this 
paper. Since all vector potential amplitudes A considered 
in this work are < 10 a.u. the dipole approximation is ap- 
plicable, as has been shown in pel]. 
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FIG. 1: The probability density, integrated over a sphere of 
radius Rt, — 60 for laser pulses of frequency u = 1 [plot (a)] 
and u = 2 [plot (b)] vs. time. (4,16,4)- and (4,32,4)-pulses 
were chosen for u = 1 and ui = 2, respectively. The E values 
are attached to the curves. The survival probability increases 
with increasing peak field strengths, but not beyond a certain 
maximum. In some of the cases an ionization rate can be 
determined from the (in the logarithmic plot) linear slope. 



A. Survival probability vs. time 

In order to calculate the probability for the hydrogen 
atom to "survive" non-ionized till the end of the pulse 
the probability density is integrated over a spherical vol- 
ume of radius Rb = 60 atomic units. The radius Rb was 
chosen several times greater than the biggest excursion 
amplitude a — 10 a.u. considered in this paper. Since the 
free part of the wavefunction vanishes inside that sphere 
at sufficiently large times after the laser pulse has passed 
by, the amount of probability density inside the sphere 
leads then with high accuracy to the same probability one 
would obtain by projecting on all bound states. The re- 
sults are presented in Fig. [l]. In the upper panel (a) the 
results for u> = 1 are presented, in the lower panel (b) 
those for u> = 2. All pulses were up- and down-ramped 
over 4 cycles. The constant part of the pulse lasted 16 cy- 
cles for u> = 1 and 32 cycles for uj = 2. By Fig. |l| the very 
existence of dynamic stabilization is proven: the survival 
probability at the end of the pulse increases with increas- 
ing peak field strength (attached to the curves) up to a 
maximum value. Increasing the peak field strength fur- 
ther yields a decreasing survival probability (not shown 
in Fig. nl). The highest survival probability in plot (b) 
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amounts to w 47%. It should be noted, however, that the 
survival probability is rather sensitive to the laser pulse 
shape and duration. This was also found in the case of 
linear polarization and is discussed at length in Ref. jio) . 

Whenever the slope of the probability curves vs. time 
is linear on the logarithmic scale an ionization rate can be 
determined. While this is straightforward for the curves 
in the lj = 2 panel it turns out to be more difficult for the 
ijj = 1-results in (a). Longer pulses of constant intensity 
are then needed to obtain unambiguous rates. 



B. Snapshots of the probability density 

For sufficiently large excursions a and adiabatic turn- 
ons one expects the system to occupy essentially the 
ground state in the time-averaged Kramers-Henneberger 
(KH) potential. This corresponds to lowest-order high- 
frequency Floquet theory [Q. The KH potential is the 
ionic potential a free electron, just oscillating in the laser 
field, would "see", i.e., Vkh(p) = V(r + ct(t)) where a(i) 
is the trajectory of a free electron in the field 






t = 59.V2 » 






x (a.u.) 

FIG. 2: Series of contour plots of the probability density 
in the rry-plane at z = during the 9th cycle of the (2,8,2)- 
laser pulse for u)i = 1 and E\ — 10. The color-coding of 
the probability density is linear. The overlayed black circle 
of radius a = 10 indicates the corresponding free electron 
trajectory. A clockwise "hula-hoop" -like Jll| dynamics of the 
probability density ring is observed. 



a(t) = f dt'A(t'). (5) 
Jo 

In the case of a single-color, circularly polarized laser 
the KH potential, averaged over one laser cycle, has az- 
imuthal symmetry with respect to the laser propagation 
direction (the z-axis in our case) and, for sufficiently big 
a, its minimum along a circle of approximate radius a 
in a plane perpendicular to it (the xy-plane in our case). 
Consequently, the KH ground state probability density 
has its maximum along that circle. By transforming back 
to the lab frame, where the ionic potential is assumed 
centered at the origin, this ring-shaped, stabilized proba- 
bility density orbits with its center along the free electron 
trajectory which itself is a circle of radius a. What re- 
sults is a "hula-hoop" -like dynamics |ll|] of a probability 
density ring about the ion which assumes the role of the 
"hula- hoop" -dancer's hip. 

In Fig. H we present snaphots of the probability den- 
sity in the xy-plane at z = during the 9th laser cycle 
of the (2,8,2)-pulse with ui = 1. The peak field strength 
was E — 10, yielding an excursion radius a — 10. Al- 
though a probability density ring, describing the previ- 
ously explained "hula- hoop" -dynamics, is clearly visible, 
the probability density is not uniformly distributed along 
that ring but has rather a sickle-like shape. Hence, the 
system is obviously not in the KH ground state but in a 
superposition of states with different azimuthal quantum 
numbers. This leads to Rabi floppings on time scales that 
are large compared to the laser period. Note that there is 
always a probability density maximum at the position of 
the ion. For lower values of a we observed more complex 
situations where no clear ring structure but two or more 
wave packets swirl about the ion. 



C. Ionization rates 

For the single-color case, accurate ionization rates were 
numerically obtained many years ago by means of non- 
perturbative Floquet theory |L3]]. Those rates are accu- 
rate in the sense that if the system somehow manages to 
occupy a single Floquet state it will decay precisely with 
the corresponding Floquet rate. Usually an adiabatic 
turn-on is needed to transfer the system from the field- 
free groundstate to a single Floquet state. The problem 
is that an adiabatic ramping of the laser pulse counter- 
acts stabilization because the system has to spend too 
much time in the so-called "death valley" of medium 
laser intensities where ionization is strong, and thus noth- 
ing would survive to be stabilized. Hence, short turn-on 
times are desirable. However, after such rapid turn-ons 
the system might be "shaken-up" into a superposition of 
many Floquet states. Under such conditions only the full 
solution of the TDSE yields reliable ionization probabil- 
ities and, if existing, ionization rates. 

In Fig. |^ our ionization rates for u> = 1 and lu = 2 (de- 
termined from plots like those in Fig. [I]) are compared 
with Floquet results |l3| and the analytical results de- 
rived by Pont and Gavrila fi m. The Pont-Gavrila rate 
for circular polarization 



was derived from high-frequency Floquet theory com- 
bined with the first-order Born approximation. In or- 
der to meet the plane wave criterion of the Born ap- 
proximation the excursion amplitude a must not be too 
small. For high frequencies and ci > 1 the Pont-Gavrila 
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FIG. 4: ATI photoelectron spectra for u) = 1 (left panel) 
and w = 2 (right panel) after (4,8,4)-pulses (for ui — 1) and 
(8,16,8)-pulses (for u> = 2). For clarity, the spectra are shifted 
vertically with respect to each other. The peak field strengths 
were (left panel, from top to bottom) E = 0.3, 0.5, 1, 2, 5 and 
(right panel, from top to bottom) E — 0.6, 2, 10 atomic units. 
For certain laser intensities excited states come into play, lead- 
ing to pronounced modulations of the ATI peaks. With in- 
creasing laser intensity the ATI peaks move towards higher 
energies. 



FIG. 3: Ionization rates vs. excursion amplitude a for (a) 
uj = 1 and (b) u> = 2 a.u. Our numerical results (+), the full 



Floquet rates Tpi from 
by Pont and Gavrila 



(solid), and the analytical result 
TpQ (dashed-dotted) are shown. 



rate should merge with the full Floquet result. How- 
ever, in |l3| ] Floquet rates were calculated only for mod- 
erate values of a [note, that, by definition, the excursion 
amplitude a for circular polarization in |l3| differs by 
a factor \[2 with Pont's and Gavrila's, and ours, i.e., 
TpQ = rpi/V2]- For ui = 2 the Pont-Gavrila rate comes 
close to the full Floquet result already for a = 1. For 
increasing a the full Floquet result oscillates around the 
Gavrila rate. In the lower frequency case u) = 1 the 
Gavrila rate overestimates ionization. Our numerical 
rates were obtained from simulations of (4,16,4)-pulscs 
for u> = 1 and (4,32,4)-pulses for w = 2. The results 
for the highest a values shown were difficult to deter- 
mine because the ionization probability curves vs. time 
displayed no clear exponential decrease. In order to de- 
termine unambiguous ionization rates also for higher a 
a more adiabatic up-ramp of the laser field and a longer 
pulse duration would be required. However, since it is 
not the goal of this paper to determine rates which are 
probably easier to obtain via Floquet theory, we do not 
proceed further in this direction. Our numerical results 
agree well with the full Floquet rates, where available. 
For uj = 1 and a > 3.5 the Gavrila rate agrees quite well 
with our numerical data whereas for uj = 2 it underesti- 
mates ionization. 



D. Above-threshold ionization (ATI) 

In Fig. H photoelectron spectra after (4,8,4)-pulses (for 
uj = 1) and (8,16,8)-pulses (for uj = 2) of different peak 
field strength are presented. They were calculated from 
the wave function immediately after thepulse with the 
window-operator method proposed in j26fl. The three 
lowest-order ATI peaks n = 1,2,3, separated by tko, 
are clearly visible although they broaden with increas- 
ing laser intensity. From lowest order perturbation the- 
ory (LOPT) an exponential decrease of the peak heights 
~ E 2n with increasing n is expected. Within high- 
frequency Floquet theory, for a> 1, (£' + nuj)a >• 1 the 
contributions of higher orders n do not vanish exponen- 
tially but only ~ n~ 2 By virtue of the first two peaks 
in the highest-£? spectra of Fig. || this trend of increasing 
importance of higher order peaks is confirmed [note that 
in the numerically obtained wavefunction after the pulse 
the ATI peak n — 3 might be slightly suppressed because 
of the absorbing boundaries] . With increasing population 
of excited states the left shoulder of the ATI peaks be- 
comes more pronounced. For certain laser intensities a 
substructure of smaller peaks appears. Since these peaks 
are separated by the energy difference of n — 2 and n = 3- 
states, that is A£ = (4" 1 - 9" 1 )/2 w 0.07, this substruc- 
ture is likely connected to the population of these states. 
However, we did not perform yet a detailed study of this 
substructure. Here we want to focus on the ATI peak 
positions. As can be easily inferred from Fig. |] the ATI 
peaks move with increasing laser intensity to higher pho- 
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toelectron energies. This is contrary to what happens 
for low frequencies and high intensities. In general, the 
positions of the ATI peaks follow 

S^=^-A cont +nu. (7) 

Here, £' Q = £q + A is the field-dressed ground state energy 
which is ac Stark-shifted by A compared to the unper- 
turbed ground state energy £q (= —0.5 for H), and A cont 
is the energy up-shift of the continuum threshold, n is 

the number of absorbed photons > n m ; n where n m i n is 

(n ■ ) 

the smallest integer that gives 

£^T > 0. For low fre- 
quencies one has A cont — U p = f-^ji where U p is the 
ponderomotive potential, i.e., the time-averaged quiver 
energy of a free electron in the laser pulse of peak field 
strength E. The pre-factor / is (for our definition of E) 1 
for circular and 1/2 for linear polarization. Since in the 
case of long wavelengths the up-shift of the continuum 
threshold by U p is bigger than the shift A of the ground 
state, the ATI peaks move towards lower energies with 
increasing laser intensities. The "drowning" of a certain 
peak no. n below £ = is the celebrated "channel clos- 
ing" in ATI physics. 

We shall see that in the stabilization regime, where 
u> > \£q\ must hold, the shift of ATI peaks is dominated 
by the ground state up-shift A so that 

(8) 

Accurate values of the field dressed energies £' in the 
stabilization regime can be obtained from Floquet the- 
ory Solving the set of coupled Floquet equations 
yields complex energies the real part of which is the field- 
dressed energy. The complex part is half the ionization 
rate, £o,fi = £ a — iFfi/2. From the high-frequency Flo- 
quet theory point-of-view £' is the groundstate energy 
of the time-averaged KH potential. In Ref. Q Pont 
and Gavrila gave KH ground state energies for H in a 
circularly polarized laser field. Because \£' \ < \£q\ the 
ATI peaks move to higher energies with increasing laser 
intensity — contrary to the high-intensity, low-frequency- 
case (fjj). 

In Fig. ^| our numerically determined positions for the 
n = 1 ATI peak are compared with Eq. where the value 
for £q has been determined (i) by the full Floquet solution 
|l3f , (ii) as the time-averaged KH potential ground state 
energy 0], and (iii) through Pont's approximate formula 

= -(In a + 2.654284)/(2?ra) J|. The agreement may 
be considered acceptable although not excellent. The 
differences are likely due to the population of more than 
a single Floquet state. Our numerical result for lo = 
1 and a = 5 lies above the analytical curves because 
the substructure in the ATI spectrum (cf. Fig. ||) also 
alters the absolute maximum of an ATI peak so that its 
determination becomes ambiguous. 
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FIG. 5: The energy E^ 1 ' of ATI peak n = 1 vs. the excur- 
sion amplitude a = E/uj 2 for u) = 1 (left plot) and lu — 2 
(right plot). The + symbols are our numerical results. The 
curves give the expected position of ATI peak n = 1 according 
(B) with £q from the full Floquet calculation |L3| (solid line), 
with £q as the ground state of the time-averaged KH poten- 
tial § (dashed), and with £' Q = -Qua + 2.654284)/(27ra) §] 
(dashed-dotted) . 



IV. TWO-COLOR CASE WITH OPPOSITE 
CIRCULAR POLARIZATIONS 

A. Time-averaged KH potentials 

It is useful for the analysis of our two-color results to 
get familiar with the free electron trajectories in such 
fields and the corresponding time-averaged KH poten- 
tials. In the two-color case the KH potentials were aver- 
aged over the low-frequency laser period 2n/uii. Let us 
consider the two laser frequencies uj± — 1 and L02 = 2. If 
the second field is absent the free electron trajectory ot(t) 
describes a circle in the xy-plane of radius cio = E^ ju)\. 
If the second field amplitude E^ equals E^\ i.e., both 
laser fields have the same intensity, the trajectory re- 
sembles an equilateral triangle with rounded edges, and 
max|a| = 5do/4. If, instead, the vector potential am- 
plitudes of both fields are equal, E^/lui — E^ /uj 2 =, 
the electron follows a convex triangle with sharp edges, 
and max|d| = Sdo/2. Finally, equal excursion ampli- 
tudes, i.e., E^/col = E^/lj2, yields a rosette-like tra- 
jectory, and max|d| = 2do. The time-averaged KH po- 
tentials look accordingly. They are presented in Fig. ||. 
In the two-color cases (b)-(d) the three-fold symmetry 
leads to three potential minima near the turning points. 
If the stabilized electron occupies the KH ground state 
the probability density will show "trichotomy" instead 
of the well-known "dichotomy" in the case of a single- 
color linearly polarized laser where two probability peaks, 
separated by 2a were observed In general, when 
ll>2 = nuji the trajectories display (n + l)-fold symme- 
try. Provided the peak field strengths are properly cho- 
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FIG. 6: Time-averaged KH potentials in the xy-p\ane at 
z = for w\ = 1, w 2 = 2. The peak field strengths E m , £ (2) 
are chosen in such a way that (a) the 0J2-field is absent, E^ 1 ' — 
6; (b) the two laser intensities are equal, B' 1 ' = = 6; 
(c) the two vector potential amplitudes are equal, E^ 1 ' /u)\ = 
E (2) /u2 = 6; (d) the two excursion amplitudes are equal, 
E {1) fw\ = E (2) /u>l = 4. In the two-color cases (b)-(d) the 
three-fold symmetry leads to three minima in the potential. 
The corresponding KH ground state probability densities will 
therefore display "trichotomy," i.e., three spatially separated 
maxima of the probability density. 



sen, for n = 3,4,... the n + 1 minima of the time- 
averaged KH potentials will lead to "quatrochotomy," 
"pentachotomy," etc. The electron dynamics in the lab 
frame should then resemble "hula-hoop" with a triangle, 
a square, a pentagon, and so forth. We would like to 
remark that the discrete symmetry of the trajectory and 
the time-averaged KH potential is closely related to selec- 
tion rules for harmonic generation in two-color fields (see 
p7[ , and references therein). If the two fields with fre- 
quencies tui and ll>2 = nu\ are equally polarized instead of 
oppositely, the KH potentials have (n— l)-fold symmetry. 
Hence, to obtain "trichotomy" uj2 — 4wi has to be chosen 
in this case. If the two frequencies are not commensu- 
rable the free electron trajectories are not closed, and the 
two-color, time-average KH potential is ill-defined. How- 
ever, we observed stabilization also in this case although 
no distinct probability density peaks were obtained. 



B. Snapshots of the probability density 

In Fig. snapshots of the trichotomous "hula- hoop" - 
dynamics for u>\ = 1, u>2 — 2 during one laser cycle (with 
respect to the u>i = 1-field) is presented. The two laser 
intensities were chosen equal so that the free-electron tra- 
jectory is a triangle with rounded edges. This triangle 
is drawn in black in each of the plots of Fig. [7| As ex- 
pected, the center of the trichotomous probability density 
(aligned along the time-averaged KH potential) moves 



FIG. 7: Series of contour plots of the probability density in 
the xy-plane at z — during the 9th (low frequency) cycle 
of the laser pulse for uii = 1 [(2,8,2)-pulse], uj 2 = 2 [(4,16,4)- 
pulse], and E\ = E2 — 10 "Trichotomy" and "hula-hoop" - 
dynamics is observed, i.e., three peaks of probability density, 
located in the corners of an equilateral triangle whose center 
moves along the black triangle representing the orbit of a 
free electron in the same laser field. The color-coding of the 
probability density is linear. 



clockwise along the black triangle. Note that the sta- 
bilized structure moves as a whole and does not rotate 
about its center. In the comoving frame, i.e., the KH 
frame, the ion swirls along the triangle, forming another 
probability density maximum. 

In Fig. y three other examples for stabilized probability 
density are presented. In panel (a), where the excursion 
amplitudes for both fields were chosen equal, no three 
distinct peaks are visible because a = 5 is too small. 
In (b) a clear example of equal intensity quatrochotomy 
(i.e., L02 = 3) is shown and in (c) a case of pentachotomy 
(ui2 = 4) where both fields had equal vector potential 
amplitudes. The density of the n + 1 peaks varies during 
the course of the laser pulse. This implies that several 
Floquet states are involved. 



C. Survival probability vs. time 

So far the two-color survival probability was not ad- 
dressed at all. If, for instance, is chosen so that 
the U2 — 2-field alone yields optimal stabilization, and a 
U3\ = 1-field with E^> in the "death valley" is added — 
what is the survival probability? In the worst case ion- 
ization would be dominated by the "death valley" -field, 
making the two-color stabilization not very attractive. 
Fortunately, it is not the field which would lead to high- 
est ionization, if applied alone, that dominates two-color 
dynamic stabilization. In Fig. |9| the temporal evolution of 
the survival probability in the two-color case is compared 
with the two corresponding single-color results. From 
panel (a) it is clearly seen that if the two laser inten- 
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FIG. 8: Probability densities in the xy plane at z = and t = 
56.84 during the constant part of a (2,8,2)-pulse (with respect 
to uji = 1). The second frequency u>2 was (from left to right) 
2,3, and 4. The laser intensities were chosen such that (a) 

& w = & m = 5) (b) ^(D = ^(2) = 10) and (c) Aw = A ( 2) = 

10. The overlayed black curves are the corresponding free 
electron trajectories along which the center of the stabilized 
structure moves. 



sities are tuned to yield equal excursion amplitudes a 
the two-color survival probability is closer to the single- 
color high-frequency result with relatively low ionization 
probability. The low-frequency field alone has a higher 
ionization probability. In the equal intensity plot (b) it is 
the other way round: the two-color result follows closely 
the low-frequency curve whereas the high-frequency field 
alone yields higher ionization because E^ is still close 
to the "death valley." Since in plot (a) all curves cross 
each other it is apparent that the final survival proba- 
bility is sensitive to the pulse length. However, it should 
be stressed that, fortunately, the two-color survival prob- 
ability is not determined by the frequency which yields 
highest ionization. This means that the two-color "prob- 
ability density shaping," as discussed in the previous sub- 
section, can be performed without drastic reduction of 
the survival probability. 

In Fig. [To] the survival probabilities at the end of 
(2,8,2)- and (4,16,4)-pulses for u>i — 1 and u>2 — 2, re- 
spectively, are plotted vs. the field amplitudes and 
E^ . Because of the probability density representing free 
but slow electrons it takes a long time until the probabil- 
ity density integrated over the sphere of radius Rt = 60 
assumes a constant value. Hence, to obtain Fig. [l0| we 
have chosen to measure the survival probability as the 
maximum probability to find the electron in a smaller 
sphere of radius R s — 7.5 within a time period of 85 a.u. 
after the pulse. The maximum probability has to be 
taken because of oscillations due to excitation of the 2s 
and the 3s state. Note that this method for measur- 
ing the survival probability slightly overestimates ioniza- 
tion because probability density which represents slow 
recombining electrons, slowly spirals back to the ion and 
does not enter the small sphere within the observation 
time. However, here we are not interested in absolute 
numbers (which depend sensitively on the pulse dura- 
tion and shape anyway). From Fig. 10 one infers that 



From Fig. 

the highest survival probability is obtained, as expected, 
in the single-color high-frequency case where £?W = 0. 
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FIG. 9: The probability to survive the laser field non-ionized 
vs. time in laser cycles of the lower frequency field with u>i = 
1. The full curves are the two-color results for L02 = 2, the 
dotted curves are the results for the low-frequency field alone, 
the dashed curves are the results for the 0J2 = 2-field alone. 
In panel (a) the equal excursion amplitude case a — 2 is 
shown while in (b) the equal intensity case E — 6 is presented. 
(2,8,2)- and (4,16,4)-pulses were used for the u>i- and U2-field, 
respectively. 



lization probability up to a maximum value, followed by 
increasing ionization if the laser intensities are increased 
further. This is different from high-frequency Floquet 
theory where the ionization rate monotonically decreases 
with increasing intensity because the "shake-off" during 
the pulse turn-ons and -offs is not taken into account. 

In our results, the field strength regime where 
EM/oji,eW/u) 2 < 1 in Fig. [To 



may be identified as 
the "death valley" of high ionization probability. Of 
course, the survival probability increases towards unity 
as E^\& 2 ' — > 0. However, this has been suppressed in 
the plot because it blocked the visibility of the results 
in the stabilization regime in which we are interested. 
Given a certain E^, adding E^> is, from a stabilization 
point of view, almost always beneficial, provided that 
is not too strong so that shake-off dominates. If, 
on the other hand, we choose a certain E^ and add 
a wi-field the survival probability may be enhanced or 
diminished, depending on whether E^ lies in its single- 
color stabilization regime or near its single-color "death 
valley." The case of a rosette-like stabilized probability 



In general, the overall behavior is an increasing stabi- density where the excursion amplitudes a are equal leads 



FIG. 10: The survival probability after a (2,8,2)- and 
(4,16,4)-pulse of frequency u)% = 1 and u>2 = 2, respectively, 
vs. the peak field strengths E^- 1 ' and E^ 2 \ The data points 
indicated by vertical lines were calculated, the shaded surface 
is an interpolation between them. Increasing survival proba- 
bility for E W ,E^ -> was suppressed for the sake of better 
visibility of the results in the stabilization regime. See text 
for discussion. 



only to a relatively small reduction in the survival prob- 
ability with respect to the optimally stabilized W2-held 
alone (around E^ = 12 in Fig. [h]). In contrast, the 
equal intensity result remains close to the low-frequency 
single-color result where — 0. Note that maximum 
stabilization is achieved for modest values of a — E/u 2 
where the probability density splitting, i.e., trichotomy, 
quatrochotomy, etc., is not yet developed. 



D. ATI in circular two-color stabilization 

Finally, let us briefly discuss the two-color photoelec- 
tron spect ra. Th e single color case has been analyzed in 
subsection HID. 

In Fig. |11| the two-color spectra are compared with 
the corresponding single-color results where either the 
ujx = 1-ficld or the lo 2 = 2-field was present. The 
laser intensities were chosen not deep in the stabiliza- 
tion regime because at such high laser intensities the 
clear ATI-structure gets corrupted by peak broadening 
and modulations (cf. Fig. |]). 

The ATI peaks in the two-color spectra are located at 
higher energies than those of the corresponding single- 
color results. This is due to the fact that the ground 
state in the time-averaged KH potential experiences a 
stronger up-shift as compared to the single-color KH po- 
tentials. It would be interesting to analyze the ATI peak 
positions in the two-color case in the same way as it was 



done in Sec. [II D for single-color fields. However, to our 



FIG. 11: Comparison of two-color photoelectron spectra (full 
curve) with the corresponding single-color results where only 
the ll>i = 1-field (dotted) and only the u>2 = 2-field (broken) 
was present, respectively. The same pulse shapes as in Fig. || 
were used. The two-color result in the left panel has been 
obtained with A = E {1) /wi = E (2) /lo 2 = 2. In the right 
panel the excursion amplitudes were equal: a = E {1) /wf = 
£< 2 >M 2 = 1. 



knowledge two-color Floquet energies for high-frequency, 
circular polarization are not yet published. 

As far as the peak strengths are concerned the two- 
color spectrum is, of course, not simply the sum of the 
two single-color spectra: the first ATI peak in the single- 
color W2-spectrum is stronger than both, the second peak 
of the two-color and the single-color wi-field. This is com- 
patible with Fig. [To] from which one infers that adding to 
the W2-field (whose intensity lies close to the "death val- 
ley") the lower frequency wi-field suppresses ionization. 
In the equal- A- result the two-color spectrum is close to 
the single-color, low-frequency spectrum. Instead, the 
two-color spectrum for equal a is qualitatively different 
from both two single-color spectra. Not only modula- 
tions appear in the two-color result which are absent in 
the wi-spectrum but also the height of the peaks differ 
significantly from each other. 



V. CONCLUSION 

We have studied the dynamic stabilization of atomic 
hydrogen against ionization in single- and two-color, 
circularly polarized laser pulses by solving numerically 
the three-dimensional, time-dependent Schrodinger equa- 
tion. In the single-color case we confirmed the "hula- 
hoop" -like dynamics for sufficiently high excursions a. 
Ionization rates were compared with those from (high- 
frequency) Floquet theory, and good agreement was 
found. The positions of above-threshold ionization (ATI) 
peaks in our numerical photoelectron spectra were ana- 
lyzed. It was found that the ATI peaks move with in- 
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creasing laser intensity towards higher energies, in con- 
trast to what happens in intense, low-frequency fields. 
This behavior can be explained by the laser field-induced 
up-shift of the initially populated, field-free ground state. 

For two-color laser fields of opposite circular polariza- 
tion we presented stabilized probability density struc- 
tures. For laser frequencies oj\ and L02 = nuj\, n = 
2,3,... and sufficiently large excursion amplitudes (n+1) 
distinct probability density peaks were observed, in ac- 
cordance to what one expects by virtue of the two-color, 
time-averaged Kramers-Henneberger potentials. As the 
generalization of the well-known "dichotomy" in linearly 
polarized laser fields, we called this multiple probability 
density peak-splitting "trichotomy," "quatrochotomy," 
"pentachotomy" etc. 

The survival probability after two-color pulses of fre- 
quencies loi = 1 and u>2 = 2 was calculated as a function 
of the two peak field strengths. As expected, highest sta- 



bilization is achieved with the high-frequency field alone 
for a certain, optimal laser intensity. However, adding 
the second, low-frequency field, even if its intensity falls 
into the "death valley," fortunately docs not drastically 
reduce stabilization, as one might expect from a naive 
superposition point-of-view. 
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